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ABSTRACT 


The present investigations study the applications 
of various loss models for axial turbines in a strearaline 
curvature method calculating scheme. Among the loss correlat 
used are the recently modified form of ■‘-^inley-Mathieson and 
Dunham-Came loss correlation, Bal je-Binsley correlation 
and its modification and a loss correlation based on Soder- 
bergs loss model. Extensive analysis has been carried out 
for the cases of single stage, three stage and four stage 
turbines including off-design calculations. The computed 
results have been compared with the available experimental 
data for these turbines. The results indicate that the 
modified Ainley-f'fethieson and Dunham-Came loss model gives 
accurate prediction of various parameters including effici- 
encies. Bal je-Binsley loss correlation appears to be equally 
good if its , secondary’- loss correlation is replaced by more 
realistic Dunham-Came loss model. 
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CHAPTER - 1 


INTRODUCTION 

All those devices in which enei'gy is transferred either 

or from, a continuously flowing fluid by the dynamic action 

of one or more moving blade rows are classified as turbomachine, 
two 

There are/main catagories of turbomachines .■ 

i) these which absorb the power to increase the fluid pressure, 
(compressors , pumps iiDthose that produce the power by expanding 
the fluid to a lower pressure (turbinesi. Itirbines are further 
classified according to the natiire of the flov; path through 
the passage of, the rotor- When the path of the through flow 
is wholly or mainly in a plane perpendicular to the rotation 
axis in a tuibine, then it is called a radial flow turbine, and 
when the :path of the through fldv' is wholly ox' mainly parallel 
to the rotational axis, it is termed an axial flow turbine* 

VJhen both the radial and axial velocity components are present 
in significant amount, the turbine is called mixed flow 
turbine. Only axial flov; turbines are used for aircraft 
propulsion. ■ 

In order to predict the performance of a turbornachine, 
we have to analyse the flow in it* The internal flow in 
turbornachine in general is nonsteady, three dimensional and 
viscous . Due to the extreme complexity of flow, the velocity 
distribution through— out the flow field can— not be calculated 



2 . 


taking account of all these factors simultaneously. The 

• 

design of the turbine consists of three major steps, the first 
is the determination of the requirement of the flow, work 
and speed, The second step is the evolution of the velocity 
diagram which should be consistent with the desired efficiency 
and number of stages. The thix'd step is the design of the 
blading that v;ill provide the necessary flov/ angles and 
velocities required by the valocity diagx'ams. 

In eax-ly days majority of axial turbines used to be 
designed by using rather simple approach in ^‘;hich the flow 
at each of the design station v.'ss assumed to have ; 

(1) Radially constant values of total temperature and total 
pressure. 

(2) Radial component of velocity are zero. 

By using these assumptions the radial equilibrium 
equation is considerable simplified. One ot the solution of 
the simple radial equilibrium equation is the free vortex 
solution in i\?hich the axial component of velocity is constant 
and tangential component varies i.nversly with the radius. 

The principal disadvantages of free vortex approach to the 
turbine design are that : 

<1) Design method is very restrictive- 

(2) Model is basically inaccurate, turbine test results 

clearly show the significant , variation in axial velocity, 
total pressure and total temperature. So it is unrealistic 
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to assume that, there* v/ould not be significant radial variations 
of total tempersture and total pressure. 

Since the advent of modern computers, more effective 
for 

methods/solving radial equilibrium equation retaining the 

have come into vogue 

radial variation terms/ The three dim.ensional flo^v problem 
is usually sirriplified by first solving the flow problem in 
various 2-D surfaces i.e. fiov/ solution on mean hub to shroud 
stream surface v/hich is commonly called a nieiiO'iional su.rface, 
the other sux'faces being blade to blade surface and crthognal 
surfaces. The blade to blade flov; solution yields the blade 
surface velocity distribution. The following are the recent methoc% 
to solve the flow problem in blade to blade flow, some of 
which are equally suited to meridional flow solution. 

STREAM line; CURVATURE i'-lETriOD 

The stream line cxirvature method is also called as 
the stream filament method. This is a numerical method for 
solving the fluid dynamic problem. The basic model used for 
st.ream filament analysis of a turbine design consist of a 
series of streamlines which trace the path of the flow from 
known condition to the first stage .and then to the final 
design station at the exit from the last blade row. Assuming 
that flow is axisymmetric through out the entire turbine, a 
series of the streamlines can be selected at tur.cane inlet 
to define a series of annular stream tubes . The positior^ of 
the streamline is defined so that the adjacent stream surface ■ 
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contains a kncvjn fraction of total flow. Since the flow is 
assumed axis NMiime trie, the flow through the turbine can be 
represented by a meriodional section as shown in fig. (i ) 

Hov/ever for, an analysis of geometry and design point 
performance of turbine, the rrain objective is to define the 
blading recfuirement also, 'turbine is subdivided into a number 
of elements. 'These individual elements can not be considered 
in isolation for these are inter related bi.' censideration 
of radial equilibrium and overall design requieeinents which 
include satisfying continuity eepation at each, of the design 
station. Hence the principal r-equirement for the ctrca:;i 
filament analysis of a turbine design is the solution of 
radial equilibrium eqiaation. This equation must take into 
account the radial variation of entropy & enthalpy. 

In the meridional plane tho flow equation may be 
written in the fonn of a gradient of meridional velocity, 
with radius. However, there are some additional terms associated 
with radial gradient in stagnation enthalpy and entropy, 
arising from work transfer and losses within the blade row. 
Additional parabolic equations for the change of the^e 
properties along the meridional stream line must be satisfied.- 
The streaml.ine curvature method of , solution is used to solve 
%his equation for velobity iteratively . The method is appli- 
cable to transonic flow as long: as the -Kach number based on 
meridiOhal velocity is nOt close to unity* 
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(2) TIME MfeRCHING T.a;CKNIQU£ 

There are several difficulties in implementing the 
streamline curvature and the matrix method in the supersonic & 
transonic regimes . 


In the time marching technique the form of the equation 
is same throughout the flow field v^hether the flov; is subsonic, 
supersonic or transonic in mrure. The basic principle of 
time marching technique is to start vjith a guessed flow 
distribution and integroite the time dependent eqxiation of 
motion and energy forvjard with time until a steady solution 
is obtained. It is 'possible to take care of the presence of 
shock vTaves by such techniques as shock fitting and shock 
capture. 


(3) MATRIX MTTHOD 


The matrix through flow method is so called because it 
calculates the flow within or through the blade row. The 
technique involves covering the region of interest with a 
fixed regular grid (as sho'V'mi in Fig. (5A) and v/riting ■ a 
finite difference approximation at every interior grid pt. 




(l.l) 


This will results in one algebi^ic equatiori for 
every grid point in terms of stream fuDction at that and; 
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neighbouring point. The system of equations can be expressed 
in matrix form as 

= 1^1 ( 1 ^ 2 ). 

jr '*X 

where is coefficient matrix derived from replacing the 

differential operator ( ), is the vector of 

unknoi^/ai stre-am function value, and fcl is thevector of 
quantities q (x,y) froiTi equation (l.l). 

Since the R.H.S. of equation (l.l) is a function of 'V" 
and its derivative, the system of equation is nonlinear' and 
must be solved iteratively i.e., by first estimating 
computing fol , and then repeatedly solving the equation ( 1 , 1) 
for [->-] . The value for is improved/each iteration 

using the previous value of 

Since A is the function of grid shape only, it nee;d 
be computed and inverted only once. This is done by factoring 
1^3, which is a square; matrix, into triangular matrics M 
and tu) . They can then be used for successive ' iteration 
and different boundary coditions. This method offers fast 
convergence/ and stability at high flow rate and machine 
speed. 

(4) FINITE ELEMENT METHOD 

This method is very Mich similar to the matrix method. 
The main differehce is that instead of^sing finite differencing 
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technique as in the matrix method, the- finite element method 
is used in solving the governing equations. This method is 
comparaitively neu and in the process of development. Hirsch 
& VJarsee ( 28 ) have described this method in much details. 

IjCSS flO'ie.ptLj « “ 

The losses arise in a tiirbine are due to the viscosity. 
Three dimensional effect (secondarx-^ flow, tip clearance) and 
mixing (interfilament m.ixing) . These losses are called as the 
profile losses and the secondary' losses. Analytical calculation 
of profile losses is difficult. in a real machine this loss 
can-not be isolated from other sources of loss, so prediction 
of profile loss are generally based on and compared with test 
made in cascade wind tunnel. There are tv;o methods for prediction 
of the profile loss. 

(1) This re^quires only simple blading parameter such as flov; 
angle, pitch chord ratio, Reynolds number and Mach dumber. 
This method is called the correlation method. 

(2) A second method has been made ;^ossible by the advent of 
the digital computer. The loss is obtained from the 
calculation of velocity distribution and boundary' layer 
growth on the blade surface. This method is called as 
calculation method. 

The correlation method is necessary when deciding 
upon optimum stage geometry' of the turbine design and would 
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usually be used in predicting the machine efficiency. The 

calculation method can bo used in place of casca.ie test when 

optimizing the blade profile. In many situations this method 

can not be used, e.g. when detailed blade profile is not 

known or because of the time required in data prepration would 

be more. The objective of vaiuous loss definition methods is, 

to provide for a loss in total pressure as normalized by some 

suitable normalizing procedure which correlates well with 

the experim.enuel evidence. The basic purpose of the loss 

model can be described in other way also, the equation of 

flow used in the meridional flow analysis are written for an 

inviscid flow. Ihe analysis of the inviscid flow would be 

of little importance in the prediction of performance of 

turbine unless some metliod is found to take into account 

change in entropy due to viscosity. This is done by using a 

loss model. The loss model superpo.--es on the inviscid flow the 

loss 

changes in entropy (Effect of viscosity) by means of a/equation 
(in the form of stagnation pressure loss or entropy change 
or loss in kinetic energy) . 

It is evident that sophisticated numerical analysis 
of a selected turbine design req^iirernent would be of little 
value if the blade row performance data used in conjuctlon 
with the analysis were not consistent with the blading to be 
used for design. The evaluation of losses in flow is achieved 
from correlation , of experimental data. A large amount of such 
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data is available in case of axial compressors. A NA.SA special 
publication (liASA SP-36') on axial compressor has a collection 
of such data. Mazurndar ( 36 ) has reviewed the different 
correlation for axial compressor. Such a review or collection 
of data is not available in e. single source for the case of 
axial turbines. There are several loss correlation available 
in literature howevexv Some of these are rather old and may 
not be particularly suited for modem blade technology. 

Some of the well known loss correlations for the prediction 
of the performance of axial flow turbines are, efficiency 
correlation of Smith (26), Soderberg’s loss correlation ( 22) 
and Ainley Mathieson loss correlation (23). The secondary 
loss correlation of Ainley Mathieson ( 23 ) was modified by 
Dunham ' &. Came ( 24 ) . Another important correlation for 
evaluating the losses is given by Balje and Binsley ( 27) based 
on Zweifelfe criterion for blade loading and Truckenbrodts 
relation for momentum thickness . 

REVIEW OF LITEl-c'-.TURE 

Internal flow in turbomachine is extremely complex, in 
general it is three dimensional, viscous and unsteady. The 
most successful approach to the problem has been to solve the 
flow in two families of intersecting surfaces. Blade to blade 
surfaces and me riodional surface. In the analysis of , blade 
to blade plane problem, ths^ flow relative to blade is 
assumed to be steady and irrotational . Such solution are 
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then usc^d in calculating the profiles losses using calculation 
methods. Solutions obtained by conformal rrans formation for 
incompressible flow ar^. given for example by Merchant Sc COLLAR 
( 1 ) and GOST BLOW ( 2 ) but this method is limited to the 
simple cascade gosometries . Tho: singularity method, in which 
blades vSre replaced by vortex and sink-source distributions 
has been used by SCHLICITNG and SCHOLL (37 ) , ISaY and 

lARTENSEM ( 3 ) and WILKINSON ( 4, ) for incompressible-- 

flows. The singularity method has . been extended for the 
subsonic oompressible flow taking into caro of compressibility 
corrections such as by VON KARiAh ( 5 ) and TSIEN ( 6 ) , 

A scheme for solution of compressible steady flow in 

the blade to b3.adc: plane was given by 'i'HEODOR.E Kxi.TSANIS ( 7 ) . 

This method is basically applicable to subsonic flow but can 
locally 

also deal with/suporsonic flows. The blades my be fixed or 
rotating. The flow may be axial, radial or miXvjd. This is one 
of tho first applications of streamline curvature (velocity . 
gradientmethod) technique. Subsequently the streamline 
curvature method used by Bindon &; Gi^RMICHAEL (3 ) in which 
a differential equation for the velocity gradient , along the 
normal to the streamline is written in terms of the radius 
of curvature : of streamline . ihls aquation is integrated 1 
across the blade passage to give the velocity profile and 
constant of integration is deterfiined by the continuity 
equation. A ; 7 WILKIHSGN: ( 9 ) used the same technique 
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but the main difference was that, th(„y used the velocity 
gradient equation along the quhsi-orthcgnal to the streamline. 
Further the same equation is \vritten in terms of the stream 
function satisfying the continuity equation, given by SMITH 
Sc KATSA^IS ( 10 ) . 

Recent development in calculation of the flows in th.; 
blade to blade includes the tim'-- marching technique. 

The difficulty cf mixed, elliptic and hytaerbolic domain has 
b£!en overcome by use of time marching techniques, treating the 
governing time dependent cgiaatic^n for compressible flow. 

These eqaations were solved by Marsh and Meriyweather ( ii ) 
in the finite difference form using a Taylor series exapnasion 
aroxind the grid point. Oaneshyar and Glynn ( 12 ) have 

developed a more efficient time m2irching scheme based on the 
method of characterstic vjhich shows comiparable accuracy v/ith 
ithat of MARSH & MERRY weather^ neither of these scheme resort 
to the introduction of stabilizing device. GOPI-.LaKR.ISHI'IAN Sc 
BOZZOIYi. (.13 ) gave an introduction of such stablizing device. 
The ■ foregoing methods do not give any explanation concerning 
shock location or entropy rise. Many flows however, contain 
only weak shocks x^j-hich may be troatod using the i sen tropic 
relation. MCDGlbtLD (14 has shown good agres-ont: foetweon : 
isentropic time marching theory and experiment for pressure 
distributi,6n in, which the maximum local Mach Number is around 
1.43. This m-rohod is ternrid a finite area method. In an 
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applicstion of finite element technique/ AlC-Y & SCEK ( 15 ) 
have analysed of the transonic flow through a cascade of 
airfoils using the finite element method. 'Ihis method consist 
of the development of a computational grid suitable for 
complex flow streture and different boundary conditions hav:e 
been discussed. Modeling of the shock and the convergence 
characterstic of the developed scheme has also been discussed. 
Further the finite element technique for steady compressible 
flo'w with the possible extension to unsteady flow is given 
by SHOMFSCfi ( 16 ) . 

In the meriodional plane the flow equation may bo 
v/ritten in the form of the gradient of meriodional velocity 
along the radius. SMITH( 17) has derived the exact radial 
equilibrium equation in a form siiitable for the stream line 
curvature method. The equation has been derived by circum- 
ferential averaging of the flov.? variables in the momentum 
equation. NOVAK ( 18 ) has described in details the process 
of solving the meriodional flow problem using the streamline 
curvature method. He has derived the governing redial equilibrium 
equation from the radial momentum equation, or an inviscid flow 
by making the use of axysirmc trie assumption . KCVAK (18 ) 

has also discussed the problem of an non-axisymmetric 
computation which is necassaiy^ v,^en epmputations are to be 
made within a blade row where the axysy.rrmetry assumption ip 
not valid. In a more sophisticated analysis, NOVAK & HEi-^aSEY 
(19) have obtained a “nearly three dimensional” solution to the 
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flov; in a turbomachine by using the streairdine curvature 
method and iterating between the solutions obtained on the 
blade to blade and meriodional planes. I'KOSi- ( 20 ) 
has .also given a: method so as to be 3£)3.e to make th€; compu- 
tations vdthin Dlade I'ov:. This lias bean achieved by using 
the concept of a “m:-.:an scroavi surface’’ developed for matrix 
method. The flow equations have been derived on this stream 
surface and the blade forces determined as that required to 
keep thfes flow in the surface. FROST (20 ) called this 

method the 'Stream line curv^ature through flow method' , In 
an application to the g.as turbine design, CiiRTEj-., PIATT and 
•LEMHERF. ( 21 ) have used tho' stn.;amline cujcvature method 

for the evaj.uation of th aosign point perf.-irmance of axial 
turbine. The radial equilibrium equation has been simplified 
by neglecting the term involving the change in momentum in 
the meriodional direction by assuming that streamlines slope 
are small. 

There are several loss correlations available in the 
literature. They are the SODSRBSRG bOSS COai^L^TTION ( 22 ) , 
Sodtsrberg has correlated the losses on a basis of sp-oce . 
chord ratio, ■ Reynolds nurriber, Aspect ratio, thickness ratio 
and blade geometry. AINLSY. & liATlISSON < 23 ) also correlates 
the profile drag at zero incidence and shows the profile drag 
for nozzle and impulse blading at various values of outlet 
angle and various space— chord ratio. D.tl'jHAfvi & GJUiE ( 24 ) 
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gave their correl3tioj| for the secondary* losses. Kh.CKSR & 
ClCkPUtJ ( 25 ) described a mean line loss system capable of 

predicting the design point efficiency of current axial 
turbines.. They modified the secondary/ loss calculation given 
by AIN uSY/Mf-iTHIESOK -u ( 24 ) . Th y describe 

the dependence, of the secondary losses on the blade aspect 
ratio. It is assumed that aocondery losses varies as 
reciprocal of blade aspect ratio, over tho complete range of 
aspect ratio. A modi fie? ^tion tO' Soderberg's loss correlation 
was devejloned by CtRTER et. al. ( 21 ) . They Have made use 
of Smith's efficiency correlation ( '26 ) in making ttiis 
rjiodif ic^tion . 31NSLEY ( 27 ) derived a generalized 

loss correl- .tion for tgrbinaa. The effects included were 
Reynolds number, bl.ade angle, blade height, number of blades, 
trailing edge thickness, tip clearance and degre-..; of reaction. . 

SCOPE OF PRESENT INVESTIGATION ' 

■ In the present work a comparitive study of axial 
turbine loss-iTK)dels in a stream line curvature commuting scheme 
has been done. Earlier Govindan (39) had used the.AMDC loss~ 
correlation, how''ever it over predicted the losses. The modified 
form of this model has recently been given by Kacker, and 
Okapttu (25), in which several correct-ion factors have been 
introduced in the profile and secondary loss corrc-lations 
and the level of profile less h'^s also been reduced. In 
addition Balje Bins ley loss model, and its modification have 
also been used. Comparision of the loss models have been 



done on the basis of design pressure ratio as well as based 
on the design mass flow. The Carter's loss model has been 
used for single soage turbine, and for multistage turbine 
for a design pressure ratio. 

Extensive computations have been carried out for 
three test turbines using these loss models. Computations 
have also been carried out for off -design case by varying 
the mass flow rate at the design r.p.m. The specific heat 
values have been specified according to the tempera t\ire 
values at the various design stations. A fifty percent 
value of mixing coefficient (for interfilament mixing) has 
been used. The use of mixing coefficient is essentially as 
a help in the convergence as its effect on computed results 
is insignificant. The results have been plotted using the 
Graphical Terminal of DEC - 10 computer. 



CHAPTER - 2 


gOIU>IlJL;^TIOH & APPROACH TO SOLUTION 

2,1 FORMULATION 

Three dimensional inviscid pattern of flow can be 
solved by using the following equations 

1. Equation of motion 

2. Energy aquation 

3. Continuity equation 

4. Equation of state. 

By solving these equation we can determine the fluid 
properties li3<e# density enthalpy entropy and the three 
component of velocity. If we combine the three equations of 
motion with energy equation, we can obtain an equation for 
reversible adiabatic flow which states that entropy of fluid 
remains constant along any stream line. If a reversible 
adiabatic flow takes place on a prescribed stream surface, 
then there exists a.- geometrical condition relating to the 
three components of velocityvector in order that the velocity 
vector lies on the stream surface. The flow pattern is obtained 
by replacing one of the equation of motion ;with the geometrical 
condition for the flow to remain on the stream surface. The 
replaced equation of motion can be used to. evaluate the force 
necessary to keep the” flow, on the stream surface. This force, 
for an inviscid flow, lies normal to the surface and is 
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t 

contained only in the equation of motion for this direction. 

Thus for an reversible '-idinbatic flow on a prescribed stream 
surface the flow pattern can be dctentiined by the fcllov/ing 
eqnations;- 

(1) - Continuity 

(2) (a) Any two of : x equation of motion 

r equation of motion . 
s entropy 

(b) Geometrical condition 

(3) Knbrgy . 

(4) State 

For an irreversible adiabatic flow, the entropy increases 
along the streamlines and this effect is taken into account 
by means of the loss model. Since the irreversibility of 
flow is cS.used by the viscous effects . The introduction of 
a loss model superposes the viscous effecrs of the real flow 
on an otherwise inviscid formuJation. 

The follov/ing one the set of equations that are widely 
used in the analysis of the flow in turbomachine and has been 
used in the present investigations. 

(1) Continuity'; ! 

(2) (a) Equation of motions 4 r direction) or (Radial equilibriumr 

equation),,' ' ■■ y , ' y ' 

(b) entropy (loss model) 

(c) Geometrical condition ; 
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(3) energy ecjuation 

(4) equation of state 


2,2 DISCRIPTIGi'i OF EQUAilOi^irS 

The follo-'.'ing assumption have been made in the derivation 

of radial equilibrium equation 

(1) Fluid is frictionless 

(2) The Rotor is rigid and rotates with constant angular 
velocity 

(3) Flow is steady and axysymmetric 

(4) The fluid is a somi perfect gas i.e. the equo. tion of 
state is P = ^RT ^2.1^ with R constant &. the specific heat 
are dependent only upon the temperature. 


For a frame of reference rotating with constant angular 
velocity about the z axis, Nexintons 2nd law of motion 

gives for frictionless fluid. 


VP ^ o nf 

= Frr - ^ + 2 ui xw 

f Dt 


2.2 


It is convenient to use the relative cylindarical 
coordinate system/ r, Q &’i/ The , unit vector derivative's are 


bl 
: 1 

Dt 


W 


^ & :,^^u 

u .. 'it“ 


w I 
u r 


and the kinematic relation 


. . r 

Dl ) 

: Dt ■ 


=: W 


DC: ' '3 

■Djp 


The capital D operator , signifies that a fluid particle 
is being followed during the differentiation. The quation (2,2) 
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can be v/ritten in berms of its scalar components are 

f W - 

( u 


1 ^ P 

f 73 r 


b:s 


(2.5a) 


1 P 

J ^5 c 




u 


Da 


'•■J 

' u 'r 


+ 2 Q W 


(2.3b) 


1 dp 


3 X Ds 


(2.3c) 


By the use of the substitution V 


V'/u+ O r, between the 


aosolute relative velocities, the equation no. (2.3a) & 
(2.3b) can be v/ritten as 

v2 

, (2.4’) 


1 P 

1 ^r 


K 


Dr' 

r 

DS 


u 

r 


1 


S’ 


-^ ( r V ) 
r Db V 


(2.4b) 


we can also write 


W' 


Di 1 

ds 


= W 


X 


( ) 

Dx 


(2.5) 


Where Dx repr -sent increase in the x coordinate the particle 
unciergoes as it moves a distance DS in the flow direction.. 
Equation like (2. 5), can be v/ritten for any desired" direction. 
^ direction of p^articular interest hepe is the meriodional ■ 
direction defined by 


I D' = I D + I D 
. m m XX r . r 


(2^6)- 
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Equation (2.5) becomes as 


( ) 

DS 


= 'JM U 

m 


( ) 


jm 


(2.7) 


Substituting equation (2.7) in equation no. (2.4a), we get 

1 32 ilr lu ■ ( 2.8) 

S dr 


m 


_u 

r 


ill r the r we can write VI = W Sin 0 

r m 


Equation (2.8) can be rewr-ituen as 

,2 


(2.9) 


Now 


and 


1 ^ P 
3 ^ ^ 

P Sln^ 


V' 


u 


Cos<$ 


'Jl 

rn 

d4? 


2 „ 3 ^ 

ra D m 

m 


D W 


m 


D 


m 


m 


d4^ 


m 


D 


m 


.( 2 . 10 ) 

(2.11) 


1 ' 


rn 


The radius 'of curvature is choosed negative when the streamline 

pro jection . on the meriodional plane is concave upwards 

/ / 

Now the equation (2.10) can be rewritten again 

r2 


1 ^p- 


DW 


W ■■ 

— — = — + Cos d? - W 

a r r r 


m 


( 2 , 12 ) 


in terms of the absolute velocity components, equation (2.12) 
becomes as 


1 ^ P ^ _ 

5 ~ r 


2 ' ■ - '2 

\ 'v P -v 

^ + Cos4» - :Sln <^ ^ 


(2,13) 


m 


D 


m 
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equation (2.13) ig written in more convenient forms by 
replacing the static pr-essure by the stagnation pi'essure 



V2 + ^ 

u m 1 

2C To { 

P -i 


from 


we get 


2 2 

u m 


1- dPo 
Po dr 


3^ 

To' 


dTo 

dr 


- V 


DV 


u 


u 


dV‘ 


m 


dr 


dr 


+ Rio 1- 


cos4> 

_ii + _nL_ — . _ V 

r r m 

m 


Sin 


DV 


m 


V^+ 
u m 

2Cp To_ 


(2.14) 


D 


DV 


Again it is necessary tocacpress the terms 


m 


m 


in terms 


of the meriodional path line slope & curvature. By making 
the use of continuity equation, equation of state, laxvs 
of thermodynamics, it can be shown in (Appendix A). 


DV 


Sin 


m 


fv, 


ni 


(1+m| -K 1 eosa 


D 


m 


i 


Sin^ 


+ ta 



(2.15) 


m 


1- M 


m 


2,2.2 CONTINUITY SQUATIOb 

The continuity equation provides the necessary condition 
for the determination of constant of integration when integrating 
the radial equilibrium equation. The location of the stream- 
lines at each of the design station and meriodional velocity 
at each of btie streamlines are determined from the mass 
flow continuity equation. For an axysymrnetric flow passing- 
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through an axial section or an arbitrary annulus, the continuity 


aquation is as follows - 

r- 

where W.„ is total mass flov; in ky./sec., ^ is the static 

'X 

density in r, & r are hub and casing radii of the 

h c 

sections . 


Reformulating the equation no. (2.16) in te.rrr;s of 
Total- Temperature, Total pressure. Tangential &: mex'iodional 
components of velocity. Substituting for the density in (2,16). 




V 


,2 




1 


(2.17) 


The equation no. (l.l6) can be expressed as follows 


W, 


T 


“F" 5 


■Po 

To 


1 - 


2 2 
u 


+ V- -T 
u m I 

2 C T^ I 
p o 


1 

X 

V cos <» r dr 
m * 


(2:18 


Since the radius dependent variable of above equation 
will not be simple analytical function of radius, the distri- 
bution of meriodional velocity which satisfy the continuity 
equation will have to be ob-tainod using an iterative numerical 

procedure. 

2.2.3 THE LOSS EgUATIOH 

Since the total temperature & total pressure are 
choosen as the principal variable in the analysis, the 
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i rreve rs ibi 1 i ty 


terras of stagna 
coefficient is 


Y 


K 


or 



or flow is more conveniently expressed in 
tion pressure loss. The total pressure loss 


defined as fellow 


Poo 


Pc. 


Pc. 


P. 


(2.19) 


Poo 

C 



( 2 . 20 ) 


For the rotor, loss coefficient is defined .as 


= 


P6^ 


J 

Pb, 


( 2 . 21 ) 


t 

Po 


2. 


or 


PO, 


j 

PO 


2s 


1 + Y^( l-Pj ) 


( 2 . 22 ) 


Po' 


The above rels-tions are used to determine the stagnation 
pressure loss when the total pressure less coefficients are 
known from the loss correlation of experimental loss data in 
the form of a loss model. 

2.2.5 THE WORK EQUATION 

The total temperature will be obtained fix)m the 
specified work which is readily^ expressed in terms of the 
total temperature, drop. Thus 

' C A To = W , 'T . ' ‘(2.25) 

p*'* , 
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w is the work extracted along the streamline the tangential 
velocity component is obtained using the buler ■work equation 
as 


C 

P 


^ To 




V- ^ u. 
u2 2 



The equation is used in analysis when design vjork 
is specified. 


2.2.6 GEOMETRIC^.L COND I TI Oh 

At a stater exit the tangential velocity is related 
to the meriodional velocity from a specification of exit flow 
angle in velocity triangle, i .e . 


V = V 
u m 


cos tan 


P 


2.2.7 


(2.28) 


Denton ( 30 ) also described stream line curvature 

method for through flow calculation fox' transonic axial flow 
turbines. Frorrt the assumption of : axial symmetry i,t is possible 
to de find a series of meriodional stream surface. %ere are 
surface of revolution along which the fluid particle are 
assumed to move through the machine. The principle of stream 
line curvature is to write the equation of motion along the 
lines roughly perpendicular to these stream surface (i.e. 
along' quasi-orthogohal lines) in terms of the curvature of the 
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surface in the meriodional plane. Qua si -orthogonal is a 
fixed line vvhich goes frori one v/all to thoother of the channel 
confining the; flov.galloi; the mass flov/ therefore crosses it 
and this is used to satisfy the continuity equation by 
adjusting the velocity level along the quasi-orthogonal to 
give the correct mass flov.-. 

2-3 APPROACH iO THE SOLUT-uOM OF THE EQuATIOK s- 

The solution of the flov7 field in case of the axial 

tux'bine essentially consists of obtaining value's of total 

temperature, total pressiire, and two componev:rs of tiie absolute 

velocity (V and V ) i.e. Fo, To, V and V are the principal 
“‘urn urn 

variables. For the solution we have the sets of equations in 
terms of these variables. Tho' set of equations, for the 
solution, would comprise of equa-tdon (2.14). The appropriate 
equation in the set of equation (2.19) to (2.22) and equation 
(2.25) to :'2.2S., 

It is obvious that the principal equation is the 
Radial equilibrium equation (2.19). This is an ordinary, 
differential equation, nonlinear in variable V^. The rest 
of the equations will be used for the evaluation of the other 
flow variables in equfition (2.14). However the streamline 
solope and curvature in the meriodional plane, 4* l/'^m, 

still necessitate consideration of derivatives with respect 
to the axial direction, x. i.e. 
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= tan“^ ( ~) (2.29) 

and - 2 

cl 

%’ dx ^ 

+ ( dr/dx)^'‘|-/2 (2.30) 

Hence the radial egni librium equation contains the 
derivative w.r. to both r & x of the analysis were to be 
performed for an axisymmetric flow in an arbitroiry duct/ it 
eould readily be extended to consider the path of individual 
streamline in the meriodional xalane • Hov7t;ver, in a turbine 
design poinx. analysis/ it is unrealistic to assurns that 
axisymmetric form of the radial e-quilibrium equation can be 
extended beyond the interblade rovi space into the blade rov7. 
ThuS/ the boundary conditions for the, meriodional streamlines 
in the interblade rov. space are indeterminate at the trailing 
edge and leading edge planes defining this space. Only the 
boundary streamlines, at the inner and outer annulus wallS/ are 
defined by the assumption that these streamlines follow the 
contours of annulus walls. In the meriodional plane for the 
slope and curvature of the flow, it becomes nocessary to 
adopt an arbitrary/ solution to the problem. The streamline 
slope & curvature in the meriodional plane will be treated 
in one of the three ways 

(a) The slope and curvature of boundary streamlines will be 
obtained from a definition of the wall contour/ and then 
both 4* and will be assumed to be linear function 

of radius dets^rirdned from values at the Ti-=nll. 
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(b) The slope and curvature- will both be specified 
arbitrarily as a function of radius. 

(c) The use of spline fitting technique. 

The use of the spline fitting tec-iniqu.es for calculating 
the stre-araline slope -and curvature have been adopted by many 
authors, baneshywar, c-:t. al. (38 ) h-Bve riid-.i a critic-ol 

assessment of this method fox' use in streamline cum^ature 
technique . The instability v;hich may occur in this technique 
is due to the oscillation of the second derivc.ti've of the 
spline curve which is used to calculate th-u streamline cunmture . 
■To overcome this problem, I>u-jESHY'VT.R et.al. ( 38 ) have 
suggested the use of two spline curves, tho first one fitted 
to points through which the streamline passes, from which 
the streamline slopes are calculated, the second one fitted 
to the calculated slope -at given axial, stations and is used 
to calculate the second rSarivutive and hence the streamline 
curvature. 

DENTON ( 30) described one more method to obtain 
stream surface slope and curvature on the basis of experience 
with earlier programs. Wliich showed that overall accuracy 
was determined more by, the empirical inputs to the program 
than by numerical scheme. It was decid-ad to use a simple 
parabolic curve fit to obtain stream surface slope and 
curvature. The axis of parabola is chosen to be- perpendicular to 
the line of joining its outer two points so that calculation 
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can be continued round ducts with up to 180° turnings. The 
parabola is a good means of obtaining the cxirvature being 
more accurate than many raore conplex eum/e fits , In the present 
analysis the spline fit tecnnique has bean used, which leads 
to the solution of a diagonal matrix, which is easily done 
on a computer. 


A numerical sclution of the radial equilibrium equation 
is an iterative process, '■‘■he continuity equation will px'ovide 
the boundary conditions for the evaluation of the constant 
of integration in the radial ecfiuilibrium equation. 


2, A THE ilFFERSMTIAL EQUATIOhb 


Since the radial equilibrium equation has to be solved 

£ ime 1 tan ecus ly v/ith two basic stream lines equations (The 

tangential-momentum equation & total pressure loss equations). 

2 

rrn - .-i I- -IT dV m 1 dPo 

The principal analysis variable considered are - -g -- - , dF~' 

dV„ 

and gg . Therefore radial equilibrium equation can be 


written as 
C 


dV‘ 


SL. + C — 

11 dr ^ 12 Po 


dPo ^ 
dr ^ 13 


dv 


u 


dr 


r-=Cl4 


(2.31) 


Where the coefficient ’^13 assigned 

values at each point in a particular design plane ance a 

value of meriodional velocity is selected, '-^‘he variables 

To, 4** i ^ constants & C are assumed to have 

P 


been evaluated 
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The differentiation of apspropriate equation (2.20) 
or ( 2 . 22 )/ (depending on the station and option being considered) 


leads - 

to a differential 

equa tion 

Of the 

form of 



2 

dV"- 

r 4. 

1 dPc 

Po dr 

+ . 
z, J 

dV 

H ^ 

dr ^24 — 

(2.32) 


K^ere coefficient C , , C C can be evaluated 

2.1. 2 2 2 2 2^~x 

iii terms of input specifications and meriodional velocity. 

Similarly differentiation of either of equation (2.27) or 
(2.28) for V_ result in a differential equation of the form of 

"U. 

c ^ dPo ( 2 , 33 ) 

31 dr 32 Po "dr ' ''33 dr “ ""34 

In this case coefficient is always zero. The terms 

involved in the coefficients .are present in Appendix (C) after 
differention of app.ropriate equation. 

2-5 TECHIE IQUE FOR SOLUTION 

The rewritting of the equations described above in the 
form (2.31), (2.32) and (2.33) helps in developing an unique 

technique for the solution irrespective of the design station 
being considerd, only the coefficient of differential equations 
differ from station to station. Wow the problem is to obtain 
the meriodional velocity distribution which simeltuneously 
satisfy the radial e^quilibrium and continuity equation. By 
choosing the initial value of the meridional velocity at one 
streamline pjosition, the local values of the coefficients of 
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the set of equations (2.31), (2.32), (2. 33) can be obtained. 

These equations are now solved for the dex'ivative dv“m/dr. 

The value of mei-iodionalvelocity at any adjacent streamline 
can be obtain6?d by using the standard Runge-Kutta-Gill method 
for solving the ordinary/ differential equation given in (Appendix-r-D^ 
Using the nev/ value of the raeriodional velocity values of total 
pressure, tangential velocity and coefficients can be obtained 
at the new streamline. Thus the derivative dV'"in/dr can be 
obtained at this new screanuine also. The process is repeated 
until the meribdional velocity is obtained at each of the 
streamlines used in the analysis. U'sirig the continuity er-iuation" 
the mass flow for the meriodional velocity distribution obtained 
is evaluated. Since the solution was obtained from an assumed 
value of meriodional velocity at one point in the flow field, 
the mass flow computed will in general differ from disigned 
specified val'ue. Hence the assimaed value of the meriodional 
velocity will have to be modified iteratively' until the 
starting value is consistent with continuity 'requirement. 

This constitute the inner loop of iteration procedure. 

The evaluation of the streamline slope and curvature 
constitute the outer loop of iteration and this can be evaluated 
by one of the method as suggested above in the section (2.3) (a) 
and (b) . In the case of 2.3(a) the slope & curvature are 
function of only radius at given design station. This can be 
completed as follows : Since initially the flow distribution 
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unknown, the initial streamline position are estimated from 
equal area foraech stream tube. Hence, sti~eamline positions 
have to be relocated after each solution of the radial equili 
brium and continuity equation until a converged solution for 
streamline position has been obtained. This completes the 
solution at a design station and analysis can be proceed to 
next design station downstream. 



CHAPTER - 3 


LOSS - MODELS 

3.1. IKTRODUCTIOH 

The equation of flovj used in the streani line curvature- 
method, are for an invisicid flow. However, the analysis of 
the inviscid flow would be of little use in prediction of 
the performance of the turbine unless some mathcd is found 
to take into account the change in exitropy d’ue to viscosity. 

This is the basic purpose of loss-model. Hence, an essential 
part of the development of a design analysis cornpuririg system 
is the development of a loss correlation which would be an 
integral part of the computer orogramme. 

Loss data is limited in case of axial turbines as compared 

to the axial compressors. But there are few most commonly used 

as 

correlations suOh/Soderberg ’ s Correlation and Ainley Kathieson 
loss correlation. These correlations have been developed using 
data obtained from the early gas turbines and steam turbines. 

Four loss models hav^e been used in this v/ork. One of 
these have been developed by Carter et. al. ( 21 ) .Second 

loss model has been proposed by Kacker and Okapuu ( 25 ) which 
is a modified form of AMDC (Ainley-Mathieson & Dunham-Came) 
loss correlation. The third loss model is due to Balje and 
liinsley (27 ) which has one of the best correlation, for profile 
losses. Fourth loss model is a combination of Balje-Binsley (27) 
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and Dunham-Came (24) loss model .in which the profile loss 
is the same as in the Balje loss model butthe secondary 
losses have been introduced as given by Dunham, and Came, All 
these correlations have been discussed in this chapter. 

3.2 DESCRIPTION OF LOSS I^ODELS 

3.2.1. LOSS MODEL _ 1 

The loss correlation was incorporated by Carter et. 
al. (21) in their turbine design programme. In the analysis 
of the blade element performance, total pressure loss 
coefficient, kinetic energ^^ loss coefficient or blade efficiency 
are much useful in terms of v/hich the row performance may 
be expressed. Evidently such a model is only used when data 
is available from experimental investigation. Use of such 
model is x'estrictive to preliminary design analysis only. 
Basically this loss model is the modification of the Soderberg 
(22) loss model. 

One of the widely used correlation (due to Soderberg) 
assumes that loss is principally a function of deflection. > / 

Hence one of the first correlation attempted was total 
pressure loss coefficient verses deflection .Certain curves 
were plotted. 
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points having the same exit angle X'^ere joined by a line; and 
value of the row reaction v/ere noted for each points The data 
are correlated to a reasonablo degree of accuracy by the 
following expression; 

Y = (0.000.7d + 0.000005 43-9 ( ) +10 (3.1) 

ex ex 


The above correlation implies that when the velocity 

V. 

ratio , falls below .45^ The level of the looses for a 

ex 

selected deflection begins to increase. This result is unacce- 
ptable, Since it is to be expected that losses will decrease 
smoothly as the overall rov7 acceleration increase at constant 
deflection. The high reaction blading with low values of 
V. /V will have the high values of exit angles. Hence the 
possibility that increase in the loss was due to the trailing 
edge blockage effect. In order to avoid the above contradiction. 
Carter et. al ( 21 ) found that use of deflection as a correlating 
parameter is clearly the reason for thc:oretically unsound 

dependence on velocity ratio given in the equation (3.1) . 

. of 

So it is logical to consider tangents/ inlet & exit angles 
rather than the actual angles in the correlation. The selected' 
correlation is of the foxin 



{tan f - tan p 

'h + 85 cos V 



ex 



if 


V 

-ja. 

V 



(3-2)<a) 
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Y = 


tan fi ^ - tan p 

(a. + a_ cos S ) 
4 0 r ex 


\ r V ^ 

S’ r = 4- = ( ^n >, 

^a, + ) j 


V 

if y r ^— a 
ex 


(3.2) (b) 


VJhere a 


f 


dependent on the 


Sg = 0.03, Sj = 
are 


to ag are the coefficients 
function using a^ = 0.055 
0.157255, and a^, = 3.6. 'Cne 

O 


. The velocity ratio 
, a = 0.15, a = 0.6, 

iL ^ 

values of these 


cons tan t/obtained from Smith correlation ( 26 } .'rhis particular 
type of the correlation is Suitable for design analysis. 


The loss model discussed above has some disadvantages 
which are given bas . elow : 

(1) The use of this loss model is limited by the availability 
of data to use th^:; pressure loss coefficient. Kinetic 
energy loss coefficient or blade efficiency as input to 
the calculation scheme. 

(2) The development of this loss modo-l' shows that the 

row velocity ratio and the blade angles are the relevarit 
parameters of the correlation. As an example. , the blade 
angle at the hub of turbine is maximum and decrease towards 
the tip section of the blade. This j.s true for most of the 

turbine, so according to this correlation the losses will 

/ ' ■ - , 

be maximum at the hub and minimum at the tip region. This 

is quite contradictory to the experimental evidences, 
experimentally it has been obsem/ed that ^he losses are 
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minimum i the central region of the annulus and increase 
towards the hub and tip region. So this correlation cannot 
correctly predict the distribution of these losses across the 
annulus and can predict only the overall efficiency and average 
losses • 


3.3.2 K?iCKSR AND OK/.PUU LOSS MODEL 


AMDC loss model has been replaced by the Kacker and 
Okapuu(25ljoss Model, because A^^DC loss model over predicts the 
losses. To analyse the distribution of losses in the annulus 
of turbine, it is essential to divide the losses into the 
principal constituents, namely, the profile losses and the 


s econda ry losses. 


A start is made using the two correlation for profile loss 
coefficient (Y ), obtained from cascade data, which are for the 
nozzle type blade ( 0) and impulse type blade ( "^ 2 ^ 

of conventional profile having the thickness chord ratio of 0.20. 
The profile losses are calculated in the Ainley mathieson 
model using the following correlation'. 



(■ 


^t^c 



(3,3a) 


0.20 
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This correlation is primarily dependent on the blade angles 
the secondary and tip clearance loss data have been correlated 
by using the concept of lift and drag coefficients and can be 
given as follov/s:- 


s c 


) 


cos-p^ 


(3,3.b) 


cos 


m 


More often used correlation for secondary losses been 
given by Dunham & Came ( 24) which is of the following form 


Y 


b^/bj^ ( cos cos cos^ p^/cosj^ 

function ( (3.4) 


where is the inlet boundary layer thickness vjhich is 

difficult 'to calculate. Therefore Dunham and Came (24 ) have 

replacing the function (of by a constant. , The value of 

constant has been taken as 0.0334. 


cos 


y = 0.0334 ( c mtf 

s bj^ cos p 


•) Z 


(3.5) 




where Z is the Ainley mathieson loading parameter. 


.(-i )2 


cos 


(3.6) 


COS 


3 B 


‘ ra 


The seconds ry^ loss was distributed across the annulus 


Is 


A- r^ + 3 r + 


(3,7) 
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A,B,C are constants, and can be evaluated by the 
follov;ing three conditions 

(1) The mass-averaged secondary’- loss obtained from (3.7) is 
that evaluated using the correlation (3.5) 

(2) The secondary’- losses are zero at m;:.an stream line. 

(3) The slope of the ais tribu'cicn is zero at the mean, 
streamline . 

Thus the total loss at any stream line is 

y = Y + Y (r) . 

P s 

This kind of parabolic distribution of secoi-dary losses 
was used by Govindan ( 39). The loss model of Kacker and 
Okapuu (25) is based on the • previously mentioned Ai^-IDC corre- 
lation. The profile losses, modified by Kacker & Okapuu (25) 
has taken into account the effect of exit mach Number and 
channel acceleration respectively. The combined effect of 
■these correlation is 

K '= 1-K^(l-K. ). Therefore the profile loss coefficient 

p 2 1 . 

Yp given by AMDC (24) becomes as Y^ = 0.918 ( Y^ 4MDC (3,75) 

where K , K, and can be calculated by the graphs given 

p 1 2 

by Kacker and Okapuu (25 ) . 

They gave the dependence of secondary loss calculation 
on the blade aspect ratio. In the AMDC loss model it was 
assumed that secondary losses varies as the reciprocal of 
blade aspect ratio, over the complete range of aspect ratio. 



39 



Kacker et. al ( 25 ) noted that AMDC loss model predicted 
a more rapid increase in losses as the aspect ratio was reduced 
in the present loss system this rise in secondairy loss (with 
a reduction in aspect ratio) is loss rapid for aspect ratio 
less than two. They gave- a loss model as 


f(A.R) = fori:^h^2 

= ^^c> 2 

Y = 1.2 Y AJ'IDC K 
s s s 

where K = 1-K_ (1-K ) 

s 3 p 

and K = 1-KAI-K.) 

p 2 1 


(3.76) 


(3.8) 


The values of and can be calculated by the 

graphs given by Kackei' and Okapuu (25 ) . This loss-model 
gives the more satisfactory result than the A.M.D.C. loss 
model . ■ 


3.3.3 BALtJE-BINSLEY LOSS MODEL 
INTRODUCTION 

To analyse the distribution of losses in the annulus of 
a turbine, it is essential to separate the losses into its 
principal constituents, namely. The profile losses and the 
secondary losses. A survey and comparis ion of methods for 
predicting the profile loss of turbine had been carried 

out by J.D. DENTON (30.) . He discussed the factors affecting 
the profile losses of turbine blades and various methods of 
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predicting the loss from simple blade parameters. The methods 
are compared with cascade results and it is found that none 
of them, predicts the loss very accurately. Balje and Binsley 
method was considered the best having a standard deviation of 
0.36 in value of (Predicted loss-experimental loss), 

Balje-Binsley derived a generalized loss correlation 
& critically compared with recently published data. Effect 
included were reynolds number/ blade angle/ blade height/ 
number of blades, blade trailing edge thiclcness/ tip 
clearance and degree of reaction. These generalized loss 
relationships are for use in optimization of turbines over a 
vide range of possible operating conditions. The correlation 
express losses as a function of key geometrical parameters. 

When the expressions for all significant losses are combined, 
it becomes possible to find the particular geometry which 
maximizes efficiency for a particular set of operating conditions. 

3.3.4 LOSS RELATIONS 
CASCADE LOSSES . 

A typical cascade loss distribution is presented in 
fig. No. ( 5 )/ which shows that two frictional loss sources 

have to be distignished in typical cascade : 

(a) The profile losses, this loss in a turbine blading is due 
to the boundary layer growth on the blade surface and due 
to dissipation in the blade wake. In a real machine this 
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loss can not be isolated from other sources of loss so 
prediction of profile loss are always based on and compared 
with the tests made in a cascade tunnel. In such tests the 
profile loss is taken to be the loss in the region unaffected 
by the secondary flow near the end v;all3 . 

(b) The end wall losses are caused by the difference of pressure 
between suction and pressure side of the profile near the end 
vralls. In addition clearance losses, due to the finite gap 
betv/een rotor blades and stationary casing. 

3.3.5 PROFILE LOSSES 

The approach for the calculation of profile losses is 
to calculate the pressure and velocity distribution around the 
profile, then to calculate the boundary layer displacement 
and momentum thicknesses which are the measure of profile 
losses. To apply this procedure, definite profile shapes have 
to be assumed first. An early approach was presented by 
Zweifel ( 27 ) , who argued that the actual velocity distribution 
around the cascade could be approximated by a rectangular 
distribution. Zweifels (27 ) Analysis was not able to produce 
an acceptable interrelation between the referred boundary- 
layer momentxim thickness and cascade geometry. Losses predicted 
by this method were very high. An alternative approach used 
in this regard is to assume a different form of surface velocity 
distribution, then calculate the boundary layer momentum 
thickness based on that distribution. The surface velocity 
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has been assumed by Balje & Binsley ( 27) to vary linearly, ' 
with camber line length from the leading edge to the trailing 
edge. The Truckenbrodt 's expression for the momentum thick- 
ness can be given as 


C © -3-2 

8 ( — )■ = C n A 


X=1 



3+2 

n 


dx 


(4.0) 


n is a constant (n = 3 to 6), which depends essentially 
on Reynolds number. To apply this equation we should know the 
variation of surface velocity C with distance x. By. making 
several assumptions this may be visualized as on app-roximately 
trapezoidal free stream surface velocity distribution. Thus, 
for a constant blade height along the flow path surface 
velocity is given by 


C = C, 


‘"sin 

?2 

Sin 

fl 


+ 


X (. 


(4.1) 


Substituting equation (4.1) into (4.0) rearranging and 
integrating 



4+2 

n 


n 

-)n+l 


Cl 

( —2 ) i— . 


r 

li- 


^2 , 4+fj 
Sin ^ n I 


1- sin ^2 
Sin 


n 

n+1 

(4.2) 
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Assuming n=4 and including the terms A & C^l / ^ 
constant ‘K* the .final relation for momentum thicKness ca: 

be written as 


r Sin 


*"2 ^ 4-5 


5/1 = id 


Sin n 


0.8 


Sin pj 
~ Sin % ^ 


(4-5) 


or e/1 


K 


l-(Sin R) 


4,5 


1- (Sin R) 


.8 


(4.4) 


where the value of constant K - 0.0021 


Now 9* = 9/1 ^ 1/^V^ 


’c Sin p 2 


(4.5) 


and = H 8* 

where H is the shape factor which can be calculated as 
follows : 


H 


1 

n+1 


3A 


+ 


fs 


5A 


3n+l 


+ 


f s 


5n+l 


+ - 


(n+l!)( 2n+l l^ /3n+l) (4n+l) ^ (5n+l) (6n+i ) 


'fs 


A 


fs 



cr 


) then 


H = 


1 

n-H 

(n+1) (2n+iy 


(4.6) 


or 


H = (2n+l) 
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Now th^ profiiG loss coefficient can be given as fallows ' 


£ 

jv> 


1 . 0 ~ 


Cos^ •^{1- + Sln^ (l- 

(1« S ) 2 . 

. .~ (4.7) 

1 + 2S±n^ ( 1 -^* - ~ ( 1 - -^) 


This expression is for the incompressible flov/ - - 


3.3.6 EN D Wi-LL. LOSIES 

The endv.-all losses may be attributed to two different 
flow machanisrns, one is thtj boundary layc,= r build up along 
the wall at both tip and hub sections as flov.-' passes from 
inlet to the exit of the cascade. The other is the Svcondary or 
cross flow effect which tends to move the fluid across the 
channel from pressure surface to the suction surface of the 
adjacent blade. The follov/ing loss correlation has been 


given for minimum loss coefficient 

1 r ^4n ® lK| - 0. 

h/c ^ = K_ < 1.0388 ( oVU ' + 0.08) (l + ' > 


^ L 


Sin 


100 


+.000337 \ 10 


[ 


Sin p. 


- ' '^.08) P 


Sin ^ ° " ’I 


160 


(4.8) 


With K =0.3 for nozzle 
e 

and 

K = 0.3 + tanh (2.857 M £1^, ) for the rotor blade, 
cr 
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3.3.7 TIP CLiIARANCE LOSS 

This loss is due to the pressure difference between the 
pressure and suction side of the profile end the leakage 
flov/ over the tip of th’S blade. The einperical relation for 
this is given as 



p D 

D h 


;in/> X R 


Q 


where P =0.0696 tanh (13 S/„„) 

^ S CR 


,0.0 BALu'E-B INSL.IY VvI'I'H DUNFAH-CAMB LOoS MG^SL( B . B . D , C . ) : ~ 


Bal je-Einsley loss model underpredicts the secondary 
losses, because the selection of this now set of loss-re:lation- 
ship is to describe the axial flow turbines v/ith optimum 
solidity. This loss relationship are primarily meant for 
optimum design of turbines over a wide range of the possible 
operating conditions. So in practice, however the secondary 


losses are likely to be highe.r than predicted by Balje Binsley 

correlation. To overcome uhis limitation, present work has 

a combination of loss-models of Balje-Binsley (for profile 

losses) and Dunham-Came, (foi' secondary losses). This correlation 
found 

has been/to give good agreement X’Jith the experimental results 
of the highly loaded multistage turbines. 
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"COi'-iPUTER PROC:rRAM" 


4.1 INTRODUCTION 

The basic equations which governs the design point 
performance of axial turbine have been discussed in the previou 
chapters. These equations must be solved by numerical methods 
which can be lengthy and ti ae consuming. Therefore, it is 
of considerable importance that they be solverl using a digital 
computer. 

A program for gas turbine analysis vras developed by 
Govindan ( 39 ) . However the coniputer deck of this program was 
not available as the original program had been modified fox' 
the analysis of steam turbirics . A considerable effort was 
spent to modify this steam turbine program for the case of 
gas turbine analysis. This program is v/ritten in Fortran IV 
language and is capable of analysing both single and multispool 
units, A maximum of three spool is allowed, and each spool 
may have up to 8 stages. The absolute and relative flow fields 
are computed at the first stator inlet, at each inter blade 
row plane, and at the final rotor exit. The effect of the 
radial variation of the following quantity are taken into 
account : Jnlet conditions, streamline angle of inclination 
and curvature, loss coefficient or efficiency, whirl velocity 
and axial velocity. The effect of the coolcnt flow and 
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interfilament mixing, and a station to station variation of 
the specific heat has been included. As additional feature, 
the program allows for the internal calculation of losses 
based on the correlations. Anv- of the correlation can be used 
by changing the value of an inde^ (IPLC) in the input data. 

4.2 DESCRIFTIOiSI OF INFUi bATA 

The general specification of turbine consists of 

1 . Number of spool 

2. Gas constant of v/orking fluid 

3. Mass flow at the turbine inlet 

4. Flo’S.’; conditions at inlet of turbine (Total temperature, 
total pressure, and flow angles as a function of radius) . 

The spool specifica'tion consist of 

1. Rotative speed 

2. Rotor exit relative cingles 

Finally, The spool analysis variable consist of 

1. Humber of stages 

2. Rotor relative exit angles 

3. Specific heats and spool inlet and each blade ro-w exit 
station. 

4. Annulus geometry and axial position of each station. 

5. Mass flow and total temperature of coolent added if the 


turbine is cooled. 
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6. Streamline value of mixing coefficient for each blade 
rov7 if interfilament mixing is specified. 

7. Whirl velocity or flow angle as a function of radius 
at each stator exit. 

8. Stream line value of the power output distribution or 
the rotor relative exit angles. 

9. Stage efficiency, rotor efficiency, total pressure loss 
coefficient or kinetic energy loss coefficient, v.^hen the 
losses are not calculated from the loss correlation availabl 
in the program. 

4.3 DESCRIPTION 01' N ORi-C.L OUTr-'UT 

The output of the program consists entirely of printed 

data. This can be listed as follows : 

1. All input data 

2. Tabulated streamline values of flow parameters 

3. Tabulated streamline values of mixed and/or cooled flow 
parameters for a blade row. 

4. Tabulated streamline values of the performance parameters 
of the stator and rotor blade rows. 

5. Mass averaged performance parameter for a stage. 

6. Tabulated mass averaged performance parameters for each 
stage of spool. 

7. Mass averaged performance parameter for a spool. 

8. Mass averaged performance parameter for the turbine. 

9. In addition streamline values of the flow parameters 
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through each iteration of the computation can also be printed 
if desired. 


The most serious problem of this program is the mass 
flow convergence, at statiocaS where the choking flow secured 
and resulted in oscillation ofthe solution and consequent 
abandoning of the calculation. The reason for these oscillation 
was too large a change in the meriodional velocity from one 
iteration to the next of the inner loop. This resulted in the 
mass flovj oscillating about the choking mass flow without 
convergence. The solution of this problem was to use the damping 
techniques on the meriodional velocity 

^0) ^ y(l) + ( (4.1) 

m m ^ ml 


( 2 ) 

where V . is the value of meriodional velocity calculated 
ml 

(l) 

for the new iteration, V value for the previous iteration, 

m 

( 2 ) 

V value used in new’ iteration which is the damped value, 
m 

The value of ^ between 0.7 to 0,8 to effectively damp the 
oscillation. 


In this program use of the rotor relative exit 
angles has been made instead of the w’ork specification. 


4,4. OVER-ALL PkOGRaM LOGIC 

The computer program is composed of a main routine 
and twenty three subroutines, A list of these subroutines 
along with their function is as follow : 

1, INPUT r- To print all the input data 
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2- STRAC To calculate the slope & curvature of the hub 
Sc casing stream lines. 

2- 5>PECHT s- To determine the specific heats at constant 

pressure, specific heat ratio Sc related parameter at each 
station . 

4 . POWER To determine the- droT-> in absolute total temperature 
across each stream lines of a rotor whan v7ork done by the 
turbine is specified. 

-• 1 To determine the drop in absolute total temperature 

across each stream line of a rotor when rotor exit relative 
angles are specified. 

strip To obtain the initial estimate of the radial 

position of each stream line. 

7. STRyAL To obtain the stream line value of the items 
required for the solution of radial equilibrium equation. 

8, VMl'ITL ; To obtain the initial estimate of the meriodional 
velocity at the mean stream line. 

5* To control the logic of the calculation of the 

meriodional velocity distribution. 

RUNKUT Contains the algorithms for the Runga-Kutta Gill 
Method for solving the 1st order ordinary differential 
equation , 

To obtain a value of the derivative of the 
square of the meriodional/ydocity . with respect to radial 
position for a specified velocity & radial position. 
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^2. Vi'"-SU3 To obtain the new estimate of meriodional 

velocity at moan stream line^ 

13= REMkIN To obtain the stream line value for the quantity 
tabulated in the output '..;hich have not already been 
obtained. 

S’R'JdP To obtain the stream line; value of quantity 
for calculation at next station. If also calcuxates the 
mass averaged value which are to be printed in the output. 

15. OUTIUT To print the. result of the calculations. 

16. FLC To evaluate th.- total pressure loss coefficients 

for each stream line using the loss model. 

17. *“ evaluate the total pressure loss coefficient 
fox' each stream line using the loss model developed 
for Balj e-Binsley loss model. 

1®' *- To estimate the total pressure loss coefficient 

from the kinetic energy loss coefficient v/hen they are 
specified . 

19. SPLIjmE To evaluate stream line slope & curve -ture by 
filting spline curve through the estimated stream line 
radial position. 

20. ’ IlAPl To perform parabolic Interpolation of a 

tabulated function of one, variable if parabolic inter- 
polation is not possible, linear Inuerpoletion or 
Extrapolation of a single value is performed. 

21. SLOPE To obtain the .derivative of a tabulated function 
with respect to an indexiondent variable at each tabular 

CENTRAL diBRARY 

I f A, K^np’j r, 

A:e.No. A 
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22. 

SlMeip, Solution of a qt 

2 t of si: 

mel tan ecus linear 


■^loebric ecfuation by Gauss 

me the d. 

'..uth pivoting. 

23. 

PLC.., Modified intox'nal 

loss cc 

rrelaticn of pressure 


loss coefficient repine in g 

lehe sue 

routine PtC,., and to 


evaluate the total pressure less coefficient for each stream 
line by using mediriod form- of the id-ii'C loss model. 


An overall flow diagram for the progr-am is given 
in figure ( 4 ) . The overall control of the calculaticn 
procedure is maintained by the main routine. 

The results obtained from using th.:. above computer 
program have been plotted by using the graphic terminal. A 
separate program is used to plot these curves. By using 
computer program each data can be taken to a very good 
accuracy. Curves plotted on a graphic terminal are more 
accurate than plotting these curves manually. 
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RESULTS AND DISCUSSION 

The computer program discussed in the previous chapter 
is used to evaluate the performance of a signle stage, three 
stage, and four stage turbines having very high stage loading 
factor. In this chapter the computed results with different 
loss models have been compared with the available experimental 
data for these turbines. 

5.1 TEST CASE 1 

For the first test case a single stage turbine was 
used (Ref. 31) . The velocity triangles and the turbine flow 
path for this turbine are presented in the figure No. 6. The 
program was arun with the specified rotor relative exit angles 
and design mass flow rate for each loss model, computed 
results are plotted, for both design and off-design calculations. 

Fig. (7) shows the stage exit total pressure distribution 
with hub to tip radius ratio. It is concluded from the graphs 
that the Carter's (21) loss model and BBDC loss model are 
close to the experimental curve as compared to the Bal je~Binsley 
and Kacker Okapuu loss models. Pig. (8) presents the radial 
distribution of outlet flow angles. This again indicates that 
Carter's (21) loss model is the closest one to the experimental 
curve than the other loss model. These curve shows that the 
effect of the loss distribution is less important in case of 
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the large hub-to-tip radius ratio turbines which have small 
blade heights. Pig. 9 shows the radial distribution of the 
total to total efficiency. In this figure the efficiencies 
obtained using the four loss models are compared.Carter* s (21) 
loss model shows the maximum efficiency at the tip while the 
other loss-models give the minimum efficiency ac the hub and 
maximum efficiency around the mean line as would be expected 
in an axial turbine. Maximum efficiency at the tip is unlikely 
to occur in a turbine. The efficiency plot shows the difference 
in the loss -distributions of each loss model. Fig. 10 shows 
the comparision of the off design calculations in which the 
total to total pressure ratio is plotted against the equivalent 
weight flow, Balje-Binsely loss model, Kacker-Okapuu loss 
model and BBDC loss model shows a close agreement with the 
experimental results. The design mass flov/ is 2.507 kg/sec 
with total to total pressure ratio of 1.90. 

The results for the experimental and for all the loss 
model have been summarised in , the following table 
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Prom the Table 1 it is clear that BBDC and Kacker- 
Okapuu loss model gives result close to the experimental 
values. Carter's loss model predicts a higher total to total 
pressure ratio, and a higher work extraction. Balje-Binsley 
loss model also gives comparable results for this case, 

Denton (30) suggested that a comparision can be made by 
keeping the design pressure ratio fixed. vJhich is an important 
design parameter for turbine. On this basis also all the loss 
models give a good comparision. 

5.2 TEST CASE 2 

In the second test case a three stage (highly loaded) 
turbine was used. The design details for this turbine have 
been collected from Ref. (34). Some of the design details 
are presented in Fig. 11 and the turbine flow path has been 
shown in fig. 12. 

In Figures 13, 14 and 15, the hub and tip pressure 
distributions for different axial stations have been plotted 
for Carter's (21), Kacker-Okapuu and BBDC loss models and 
each one is compared with the experimental results separately, 
Kacker-Okapuu and BBDC loss model have very close agreement 
with the experimental results. 

Pig. 16 presents the 'total to total efficiency 
distribution with percentage blade height. This figure 
compares the : , loss -model used by GovindanC39) with the 
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loss models of the present work, Kacker-Okapuu and BBDC 
loss model show better agreement v/ith the experimental results. 
All loss models give the maximum efficiency around the mean 
line as v/ould be expected in an axial turbine. This figure 
gives an interesting information regarding the AliDC loss 
model used by Govindan (39), While the original AI-TOC loss 
model underpredicts the overall efficiency, the modified AMDC 
loss model of Kacker and Okapuu overpredicts by 2%, Fig. 17 
shows the axial velocity distribution at the entry of the 
last stage of the turbine. This distribution is parabolic. 

Pig. 18 presents a plot of equivalent specific work 
against the total to total pressure ratio. Kacker-Okapuu and 
BBDC loss models are very close to the experimental results. 

Bal je-Binsley loss model under predicts the losses, while the 
modified Bal je-Binsley loss model (BBDC) gives a good agreement. 
The carter's (21) loss model failed to give a conversed 
solution for highly loaded three stage turbine. The solution 
is sensitive to the radial distribution of losses. The results 
for this loss model have been taken at reduced mass flow rate. 
This predicts the lowest value of the specific work for a 
given total to total pressure ratio. In figure 19 the equivalent 
weight flow vs. total to total pressure ratio has been plotted. 
The trend of curve for Carters loss model (21) indicate that this 
loss model is hot suited for multistage turbines although it 
may give reasonable results for high hub to tip : radius ratio 
single stage turbines. Again the Kacker-Okapuu and BBDC loss 
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model have the closest agreements with the experimental 

results. For the same design mass flow the Kacker and Ofeapuu 

loss model gives slightly higher total to total pressure ratio. 

to 

Bal je— Binsley loss model has a lower value of total/total 
pressure ratio for the design mass flow. Table 2 compares the 
predicted overall performance parameters with experimental 
results . 


P.T.O. 
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In the table 2 , Balje-Binsley loss model indicates 
that it underestimates the loss because it is giving the 
very high efficiency lowest total to total pressure ratio and 
lowest work extraction. Kacker and Okapuu shoxirs a very close 
result to the experimental result. A comparision on the basis 
of the design pressure ratio shov/s that SEDC loss mcuel has a 
good agreement with the. experimental results. Design mass 
flow has been reduced by 0.7% in this case tc get the design 
pressure ratio. For the Carter's loss model design mass 
flow has been reduced by 6% to get the design pressure ratio. 
This model predicts very lov/ efficiency. 

5.3 Test Case 3 

In this test case^ A four stage turbine with very 
high stage loading v/as used. Design details for this turbine 
have been taken from Ref. (33). Some of the design details 
are given in the fig. (20) and fig. (21). 

The fig. (22) describes the radial variation of the 
stage exit total pressure for the design mass flow rate. This 
shows that for the same radius ratio, each loss model has 
a quite different value of pressure ratio. liirther in Fig. (23) 
same graphs is plotted with design pressure ratio. Here for the 
aame value of the radius ratio every loss model has a very 
close yalue of pressure ratio to each other. The BBDC loss 
model and Kacker-Okapuu loss models have almost the same 
value of the pressure ratio for the same radius ratio. Fig. (24) 
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gives the radial variacicn of the outlet flo\-v angles. In 
most cases, there is a good agree-.ront v'itn the experimental 
results . 

lug. 25 compares, the equivalent specific veork Vs 
total to total pressure ratio. For the design pressure ratio, 
BBDC loss model c^nd Kacker-Gkapuu loss rno'.lal gives very 
close values of the equivalent specific v.;orK;, which is almost 
matching with the c|5perimental results. Saljc-fiinsley loss 
model gives a higher value of the specific work than the 
experimental value. Carter at. al. loss model is not able 
to give the solution for desiejn mass flow rate due to 
coxivei'gence problems. 

Fig. 26 explains the comparision of various loss 
models with experimental results for off design calculations. 
This figure shows the variation of equivalent v/aight flow 
with total to total pressure ratio. In general there is a 
good agreement v/ith experimental results. For the design 
mass flov/ rate E.BDC loss model has the maximum value of total 
to total pressure ratio. 

Fig, 27 gives the variation of the efficiency v/ith 
total to total pressure ratio. This shows that for the design 
pressure ratio, present loss nK>dels have higher .valu^of 
efficiency than the experim.:.rtal values. Bal je-Einsley loss ^ 
correlation shows the maximum value of efficiency, this 
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clearly indicates that this loss model underpredicts the 

* 

losses . 

Fig. 28, Pig. 29 and Pig. 30 present stage to stage 
hub and tip static pressure distributions in the turbine - 
In most cases there is a good predicciori of tho sb-.tic 

pressure distribaition. Table 3 compares tho predret-iu 
overall performance parameter* vhth expex'iment'- i results. 

Table 3 shovjs a cemparision on the basis of the 
design pressure ratio fer this tesu case give tiiau BBpC and 
Kaker-Okapuu are the besu loss mC'I-.ls. the BrOC loss model 

there is a 1% reduction in the design mass flow, bhile in case 
of the Bal je-Binsley there is 4.5% increase in the design 
mass flow in order to get the design pressure ratio. For 
Car-ter's loss model there is 8% reduction in mass flow. 
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CONCLUoICNS 

1) Jour different loss models have been used in a streamline 
curvature computing scheme tc study their suitability in 
predicting the performance of axial turbines. The computed 
results have been compared with experimc?ntal results for 
three test cases which include three and four stage turbines. 
The results h- ve been obtained both at design mess flow rate' 
and design pressure ratio. Off-design perfornvance is also 
studied by varying the mass flow rate at the design r.p.m. 

2) The modified hMDC loss correlation by Kacker and Okapuu 
has been found to give results in good agreemeut v/ith 
experimental values. The predicted efficiencies are slightly 
(2 to 3%) higher than the oxperimental values. The or-iginai 

AMDC loss model used by Gcivindan (39) under predicted, the 
efficiencies. 

3) Bal je-Binsley correlation tends to predict low loss levels 
and therefore higher efficiencies. This is .taainly due to the 
low values of secondary losses predicted by this model. However 
v/hen Balje-Binsley profile loss correlation was combined with 
Dunham and Came secondary loss correlation (BBDC los Modal), 
the predicted results were closest to experimental results. 

4) Carter’s loss model is mainly suited for single stage 
turbines with high hub to tip radius ratio. 
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Appendix A 


An expression for the term D V in the radial ecfuilibrium 
equa tion ™ 


The continuity equation can be written as 


W D 


+ V. W 


0 


% Ds 

In cylindrical co-ordinates 
1 


V. W = 


l(r VI ) 7-W aW 

r u X 

hr rdQ ox~ 


(Al) 


(A 2) 


Also, the derivative of a quantity in the direction of flow 
(the s-direction) at any instant in time may be expressed by 
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I _ _r 

l>s W * Dr 
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(A3) 


or re-arranging. 
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X X 

Use equation (A4) in (A2),-to eliminate the x derivative. 


We get 


V. W = W, 


Since, 
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w. 
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■W:;, 


= tap and 


u 


W 


= tan? 


equation (A5) can be rewritten as 
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'5. a = \ ft • 0.6) 


Ir 


rlie 
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It is known that 

eLi = + w 

Dt ^t OS 

T . ■ D( ) _ D( I 

Along with the kinematic relation j-jg 

(A7) can be written 
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, equation 
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Equation (A8) can be used for the last term of equation (A6) 
Also since 


W = W cos4^ 
X m 


(A9) 


we get. 
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We also know that 
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Substitute equation (.12) in (A6), we get 


V - w = W. 
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\ 1 1\ ( r tan '-ft,) 
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for an axisymmetric, steady flow 


= W, 


^(r tan ) . T.T sec^ 
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(7.14) 


next apply the momentum equation in the flow direction 
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Using equations (A8) and (All), (/vlS) can be written 
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This can be rewritten as 
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For an axisymmetric, steady flow 
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(A18) 


Also, 


Dp/Ds 


(A19) 


This relation along with equation 


(Al), (A14), and (A18) can 


be used to obtain an expression tor DW^ that does not 


contain 


W, Dp/Ds or D-p/Ds. get 
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Using the definition of Mach number, M = 'W /a and M = V /a 
^ m m © u 

in equation (A20) and rearranging, we get 
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Since W = V , equation (A21) can be rev/ritten as 
m m ^ 
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Appendix B 

Evaluation of the 1:0131 ’temperature drop 
across a stage 


V'Jhen the total power output of a stage Is specified, 
the total temperature distribution at a stage exit v;ill have 
to satisfy both the specified power output and it.s distribution 
across the annulus. Since initially the distribution of mass 
flow throughout the annulus is unknov/n until the distribution 
of meridional velocity has been established, the power 
distribution is specified by nondimensional power functions 
versus the nondimensional mass flow function, w(r), defined 
below 

2 nj C r dr 


w(r) 



(A23) 


If the total power output specified is HP, (horsepower) 
the total temperature drop A through the rotor must satisfy 
the equation. 


HP^ = 5.6925 
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r ^ 
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A" 


dw 


(A 24) 


Normalizing equation (A24) with respect to the total power 
and the total mass flow, leads to a definition of the power 
function P(w(r)) expressed as 


P(w(r)) = 5.6925 


w( r) 

J 




o 


(A25) 
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Differentiating equation (A25) with respect to the non- 
dimensional mass flo'w function yields an expression for the 
total temperature drop, for the jth streamline. 


HP. 




5.69 2^ 


c 

■i' p 


dp(v;(r) ) 
dw(r) 


(A26) 


The power function versus the mass flow function , will be 
a basic spe^cification for power distribution from which the 
total temperature drops are obtained. 


When the rotor relative exit angles are specified, the 
total temperature drops are obtained from> the Euler work 
equation 


j = L- 


fv . u. - V „ 
ul 1 u2 2 


(A 27) 



71 


Appendix C 

Coetficients of the differential equations 

To simplify the logic of the computer program, a 
standard procedure has beer, adopted for the solution of the 
flow field at each inter blade row station. The different 
types of stations and the various optional specifications are 
taken into account by modifications to the twelve coefficient 
appearing in the three differential equations (Equations 2.3C 
2.31, and 2.32). This appendix presents these coefficients. 
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The total pressure would be a specified function of radius . 
Hence, the coefficients of equation (2.31) are 
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dY_, 


Again, assuming that the derivative can be expressed 
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When the work output of the turbine is specified, is 

in effect a known function of radius. Thus, 
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Vvhen the rotor relative exit angles are specified, the 
coefficients of equation (2.32) are similar to those at the 
stator exit. 
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Appendix D 

The Runge-Kutia-Gill method for the solu-tion of ordinary 

differential equations 

The differential equation is of the form 
dV^ar = f(r, V^) 

where f(r,V ) is obtained from the simultaneous solution of 
m 

the radial equilibrium equation and two subsidiary differential 
equations . 

Given the value of the meridional velocity at one 

streamline, (V .) the unknown value at the adjacent stream- 
mi j 

line, (V . , is determined in four stages : 

' mi k 

where 

0 initially 

q = f 

^ '•q^ subsequently 

= =^ik - '^ij 

<2) = ‘V'kl + ‘1- '’^2 - 'll’ 

where k^ = hf £ ^'^mi’klj 

■Sj = % + 3 [_ - -‘Jq’] - *S h 
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